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ABSTRACT 

The major morphological features of a galaxy are thought to be determined by the assem- 
(— I \ bly history and net spin of its surrounding dark halo. In the simplest scenario, disk galaxies 

form predominantly in halos with high angular momentum and quiet recent assembly his- 
tory, whereas spheroids are the slowly-rotating remnants of repeated merging events. We ex- 
plore these assumptions using one hundred systems with halo masses similar to that of the 
Milky Way, identified in a series of cosmological gasdynamical simulations: the Galaxies - 
Intergalactic Medium Calculation (GIMIC). At z = 0, the simulated galaxies exhibit a wide 
variety of morphologies, from dispersion-dominated spheroids to pure disk galaxies. Surpris- 
ingly, these morphological features are very poorly correlated with their halo properties: disks 
form in halos with high and low net spin, and mergers play a negligible role in the formation 
of spheroids, whose stars form primarily in-situ. With hindsight, this weak correlation be- 
tween halo and galaxy properties is unsurprising given that a minority of the available baryons 
(~ 40%) end up in galaxies. More important to morphology is the coherent alignment of the 
angular momentum of baryons that accrete over time to form a galaxy. Spheroids tend to form 
when the spin of newly-accreted gas is misaligned with that of the extant galaxy, leading to 
the episodic formation of stars with different kinematics that cancel out the net rotation of 
the system. Disks, on the other hand, form out of gas that flows in with similar angular mo- 
mentum to that of earlier-accreted material. Gas accretion from a hot corona thus favours disk 
formation, whereas gas that flows "cold", often along separate, misaligned filaments, favours 
the formation of spheroids. In this scenario, many spheroids consist of the superposition of 
stellar components with distinct kinematics, age, and metallicity, an arrangement that might 
survive to the present day given the paucity of major mergers. Since angular momentum is 
acquired largely at turnaround, morphology depends on the early interplay between the tidal 
field and the shape of the material destined to form a galaxy. 

Key words: Galaxy: disk - Galaxy: formation - Galaxy: kinematics and dynamics - Galaxy: 
structure 



1 INTRODUCTION 

Galaxies exhibit a wide variety of morphologies, from spheroids to 
disks to bars to peculiar galaxies of irregular shape. Many physical 
properties, such as gas content, average stellar age, and the rate of 
current star formation, are known to correlate with morphology. Of 
such properties, the one that seems most tractable from a theoret- 
ical perspective is the relative importance of organized rotation in 
the structure of a galaxy. This is commonly referred to as the disk- 
to-spheroid ratio, since stellar disks are predominantly rotationally- 



flattened structures whereas spheroids have shapes largely sup- 
ported by velocity dispersion. 

smce iHubbld ( Il926h published his original morphological 
classification scheme, our understanding of the provenance of these 
two defining features of galaxy morphology has been constantly 
evolving. Spheroids were once thought to originate in the swift 
transforma tion of an early-co l lapsing, non-rotating clou d of gas 
into stars teggen et al.l 1 19621 ; [Partridge & PeeblesI Il967l : iLarsonl 
1 1974) . whereas disks were envisioned to result from the collapse of 
clouds with high angular momentum and inefficient star formation 
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teggen et alJl962l : lLarsonl 19761) . The role of merg ers as a possible 
transformational mechanism was championed by iToomrd jl977h 
and gained momentum as the hierarchical na ture of structure (and 
hence, ga laxy) formation became accepted jWhite & Rees|[l973 : 
iFrenk et al. 1985). 

Further development of these ideas led to a broad consensus 
where disks are thought to form at the center of dark matter ba- 
les as a consequence of angula r momentum conservation during 
the d issipational collapse of gas jpall & Efstathioulll980l : lMo et al.l 
Il998i) . w hereas spheroids r esult predominantly from merger events 
(see, e.g., ICole et id]|200d and references therein). Morphology is 
thus a transient feature of the hierarchical formation of a galaxy: 
a disk galaxy may be transformed into a spheroidal one after 
a major merger, but could then re-form a disk through further 
gas accretion only to be later disrupted again by another merger. 
Early galaxy formation simulations gave a visually-compelling 
demonstration of this scenario, galvanizing support for it (see, e.g., 
[steinmetz & Navarro 2002). 

This consensus view has been broadly implemented in semi- 
analytic models of galaxy formation, where the properties of galax- 
ies are deduced directly from the physical pr operties and assem- 
bly histor y of their surrounding halos (see, e.g.. lCroton et alj2006l : 
iBower et al. 2006; Somervil le et al. 2008). For example, most mod- 
els assume that the specific angular momentum of galaxies and ha- 
los are similar, and that the merger history of the halos dictates that 
of the central galaxy. 

Recent developments, however, have led to revisiting some of 
the assumptions of the simple scenario outlined above. For exam- 
ple, it has become clear that major mergers are rare, and there- 
fore probably not the primary formation mechanism of bulges 
and ellipticals. Instead, "d i sk instabilities" tefstathioueTakll 19821 : 
IChristodoulou et alj 1 19951 : | Mo et al.l Il998t) , as well as repeated 
minor encounters, are n ow claimed to b e the main formation 
path of spheroids (e.g.. Parry et al 1 12OO9I : iHopkins et all I2OI0I : 

etali 



iDe Lucia et alj|201 ll ; lBoumaud et alj|201 ih . This has helped to al 



leviate some tension between the observed evolution of the early- 
type galaxy p opulation and the major-merg er rates predicted by 
theory cBundv et al.ll2007l : lOesch et al.ll2010l) . However, questions 
might remain open, a s the estimation of merger times from obser- 
vations is non-trivial fcotz et alj2"oTlh . 

Further scrutiny has come from direct simulation of hierar- 
chical galaxy formation. Conserving enough angular momentum 
during the hierarchical assembly of a galax y to form a realistic 
stellar disk has been challenging (see, e.g.. iNavarro et al] 1 19951 : 
iNavarro & Steinmet z' 1997), as has been pinning down the effect 
on morphology o f repeated merging, especially between gas-r ich 
galaxies (see, e.g.. [Robertson et aU 20061 : iGovemato et alj2009l) . 

The inclusion of energetic feedback, needed to prevent the for- 
mation of too many faint or overly massive galaxies, has added 
an extra level of complexity to the problem, with a number of 
studies showing that morphologies can be radically altered when 
even modest changes in the strength of feedback or its implementa- 

I • ' " Tl 

tion are mtroduced (Okamoto et al. 2005; Scannapieco et al. 2008!; 
Ceverino & Klypin 2009; Sales et al. 2010; Agertz et al. 2011; 
Brook et al .i boill : IPiontek & Steinmetzll201lh . Moreover, the den- 



sity threshold assumed for star formation also changes the coupling 
between the stellar winds and the surrounding gas, playing also 
a maj or role on the properties of simulated galaxiesl lGuedes et al.l 
|20 



naic 



( IWhite&Fr"eii3ll99lh and may be collimated by the filamen- 
tary structure of the cosmic web, espec i ally in low-mass systems 
and at high redsh i ft jKeres et al.l |2005| ; iDekel & Bimboinj|2006l : 
Ivan de Voortetal] |201lV This complex accretion geometry has 
been hypothesized to promote the formation of disks by feeding 
high angular momentum material directly to forming galaxies (see, 
e.g.. lDekel et al.ll2009l : lBrooks et alj|2009l) . 

Further theoretical progress demands increased sophistication 
in numerical and semi-analytic modeling. From the simulation per- 
spective, most studies have focussed on individual systems picked 
according to what the authors believe would facilitate the forma- 
tion of a galaxy of predetermined morphology; for example, a re- 
cent major merger to study ellipticals (e.g., Meza et al. 2 0031) or a 
quiet, rapidly-rotating halo to study spirals (e.g.. .Abadi et alj2003l : 
lOovernato et al.ll2007h . Note that this presupposes the morphology 
of the resulting galaxy, and often results in the tuning of star for- 
mation and feedback parameters until, unhelpfully, results match 
prejudice. 

Statistically-significant samples of galaxies selected in an un- 
biased way and simulated at high resolution a re needed for new 
insights, a goal that, despite valiant efforts jCroft et ^ l2009l : 
ISaleset alj|200^ 12010^. has so far proved beyond reach of even 
the fastest computers and best algorithms. The situation, how- 
ever, is starting to change, with the advent of simulations of vol- 
umes lar ge enough to include dozens of we l l-resolved ~ Li, 



galaxies jCrain et al 



20091: ISchave et alj|2010l : fHahn et al.|[201ol : 



ICen & Chisaril201 iklvogelsberger et al.ll201 ll) . 



More recently, the mode of gas accretion has been recog- 
nized as playing a potentially crucial role in galaxy morphol- 
ogy. Gas can flow to galaxies largely unimpeded by shocks 



We explore these issues here using the GIMIC gasdynam- 
ical simulation series jCrain et akl l2009f) . GIMIC targeted sev- 
eral carefully selecte d regions from the Millennium Simulation 
dSpringel et"S] |2005l) in an attempt to maximize the resolution 
of individual galaxy systems while at the same time sampling a 
cosmologically- significant volume. The first ana l yses of 2: = 
GIMI C galaxies jCrain et al.l2010l : lFont et al.l201 lllMcCarthv et all 
I2OIII) show that they are fairly realistic, so we feel confident that 
we can use them to gain insight into the origin of galaxy morphol- 
ogy- 

This paper is organized as follows. In Sec.|2]we present briefly 
the numerical method and simulations. We analyze the morpholo- 
gies of simulated galaxies and their origin in Secs.[3]and|4l respec- 
tively. We summarize our main conclusions in Sec.|5] 



2 NUMERICAL SIMULATIONS 

The " Galaxies-Intergalac tic Medium Interaction Calculation", 
(GIMIC ; rCrain et al]|2009l) . simulation series follows the evolution 
of five nearly spherical regions of radius ^ 2Qh~^ Mpc each , 
selected from the Millennium Simulation jSpringel et alj l2005h . 
These regions were selected to sample environments of differ- 
ent density, deviating by (— 2, — 1, 0, +1, +2)(t from the cosmic 
average, respectively, where a is the rms mass fluctuation on 
20/i~^Mpc scales. The regions are spherical at z — 1.5, and are 
simulat ed using the s t andar d zoom-in technique described in detail 
by, e.g.. |Power et al. We provide here a basic summary of 

the main chara cteristics of t he GI MIC project, and refer the inter- 
ested reader to ICrain et al. I ( I2OO9I) for a more comprehensive de- 
scription. 

GIMIC uses a modi fied version of G ADGET-3, a development 
of the GADGET-2 code dSpringell |2005[) that includes new mod- 
ules to treat radiative cooling, star formation, chemical enrich- 
ment, and energetic feedback. Radiative cooling is implemented 
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on an element-by-element basis and thus cooling rates evolve self- 
consistently as a functio n of redshift, gas den sity, temperature and 
chemical composition dWiersma et alj l2009h. The runs al so in- 
clude a uniform ionizing background ( Haardt & Madaul200lh . with 
hydrogen- and helium-reionization redshifts of z = 9 and z = 3.5, 
respectively. 

Cold gas with densities exceeding uh = 0.1 cm"'' becomes 
eligible for star formation and is assumed to follow an effective 
equation of state, P oc p*^'^, in order to minimize numerical ar- 
tifacts in poorly-resolved regions i Schave & Dalla Vecchia, 20081) . 
Stars are assumed to follow a Chabrier IMF ( IChabriej|2003ir and 
to form at a rate that depends on the local gas pressure and that 
matches the Kennicutt-Schmidt law (Kennicutt 1989, 1998). 

Chemical enric hment is modeled as described by 
IWiersma et al.l ( l2009h . and tracks the synthesis of 11 individ- 
ual elements. As massive stars explode as supernova (SN), they 
inject energy and metals into their surroundings. This feedback is 
implemented, in practice, by using a fraction fsN of the energy 
released by SN in order to modify the velocity of a few (rj^) 
neighboring gas particles by introducing a velocity " kick" of 
magnitude Vw to each jPalla Vecchia & Schavd l2008h . These 
parameters are set in GIMIC to fsN = 0.8, rj^ = 4 and Vn, ~ 600 
km/s, which results in a good match to the peak of the glob al star 
formation rate density jCrain et alj2009l ; ISchave et al.ll20ldb . 

All GIMIC runs adopt the same cosmological parameters as 
the original Millennium Simulation, which were chosen to be con- 
sistent with the WMAP-1 constraints: = 0.25; f^A = 0.75; 
nb = 0.045; = 1; as = 0.9; Ho = 100 /i km s"^ Mpc"^ 
h = 0.73. 

The particle mass in the simulations is 1.4 x 10^ Mq and 
6.6 X 10^ Mq for the baryons and dark matter, respectively. 
The gravitational softening is initially fixed in comoving units, but 
is fixed at 2 = 3 and thereafter to e = 0.5 h^^ kpc (Plummer 
equivalent) in physical units. We shall focus here on the two GIMIC 
regions that have been run to 2 = at this resolution: the —2a and 
Oa. As we discuss below, aside from the expected difference in the 
number of systems of given mass, we see no systematic dependence 
of our results on the overdensity of the region, which may therefore 
be thought to apply to average r egions of the Universe . 

W e have used SUBFIND dSpringel et all 1200 ll ; iDolag et al] 
I2OO9I) to identify galaxies in the high-resolution regions of the 
GIMIC runs. We shall only consider in the analysis the central 
galaxies of halos within a narrow range of virialj mass: 0.5 < 
Maoo/lO^^ft-^Mo < 1.5. This ensures homogeneity in the set 
of systems selected for analysis and eliminates complications that 
may arise from considering satellites of larger systems. At redshift 
2 = 0, these criteria identify 38 and 62 galaxies in the ~2a and 
Qa runs, respectively. Each of these halos is resolved with roughly 
200,000 particles (dark plus baryonic), allowing for a reasonable 
estimate of the relative importance of the disk and spheroidal com- 
ponents. None of the results we discuss here show significant de- 
pendence on which region we consider, so we will group the 100 
galaxies together without making distinction regarding the GIMIC 
run where they were identified. 
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Figure 1. Left: The kinematic morphology parameter, K^ot ■ defined as the 
fraction of kinetic energy in organized rotation (eq.[T), versus the fraction of 
stars with circularity parameter tj > 0.5. The cumulative fraction is shown 
with a dashed line. The shaded region (ftrot > 0.7) indicates where "disk 
dominated" galaxies lie in this plot. Right: The distribution of circularities, 
e j = jz/jclrc{E), is shown for three galaxies with diiferent values of Krot ■ 



^ Virial quantities throughout this paper are defined at the radius enclosing 
200 times the critical density for closure. 



3 SIMULATED GALAXY MORPHOLOGIES 
3.1 Morphology estimates 

As discussed in Sec.[T] we shall adopt a somewhat narrow defini- 
tion of morphology based on the importance of ordered rotation in 
the structure of a galaxy. Although we refer to this as the ratio of 
disk to spheroid, it should be noted that this may differ, at times 
substantially, from traditional bulge-to-disk decompositions based 
on photometric data. The latter are based on assumptions regard- 
ing the shape of the brightness profile of disks, usually assumed 
to be exponential, and spheroids, assumed to follow either de Vau- 
coul eurs or Sersic profiles. A s discussed bv lAbadi et al ]l l2003h (see 
also IScannapieco et al.ll201Q) . these assumptions are only weakly 
fulfilled by simulated galaxies, and kinematic decompositions can 
give rather different spheroid-to-disk ratios than photometric ones. 
Photometric studies can also be affected by color gradietits, extinc- 
tion, and projection effects (see, e.g.. iGovernato et alj|2009h . We 
avoid these complications by focusing our analysis on kinematic 
data alone, although we plan to consider the implications of our 
results for photometric studies in future work. 

The importance of ordered rotation may be clearly appreci- 
ated from the distribution of the stellar orbital circularity parame- 
ter, ej = jz/jciic{E), defined as the ratio of the component of the 
specific angular momentum perpendicular to the disk (i.e., aligned 
with the total angular momentum of the galaxy), j^, to that of a 
circular orbit with the same binding energy, jciic{E). Defined in 
this way, ej takes values in the range (—1,1), where the extreme 
values correspond to the counter- and co-rotating circular orbits in 
the symmetry plane of the galaxy, respectively. 
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Figure 2. Illustration of the structure of four galaxies in our sample with increasing degree of rotational support (left to right). The first and second rows show 
edge-on and face-on projections of the stellar distribution. The yellow circle marks the radius, r-gai = 0.15 r200, used to define the galaxy. 
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Figure 3. Top and bottom rows show the circular velocity profile, Vc{r), and the stellar surface density profiles, respectively, of galaxies A-D in Fig.|2] Red 
lines in the bottom panels indicate the stellar half-mass radius of each galaxy. Dashed magenta lines indicate de Vaucouleurs' i?^/* profile fits. Straight dashed 
lines in blue indicate exponential profile fits. The shaded area highlights a radius equal to three Plummer-equivalent gravitational softening scalelengths. Note 
that spheroids can be well fit by a single i?^/^ law, but that disk-dominated galaxies show evidence of a de Vaucouleurs' spheroid plus an exponential disk. 
Disks are extended and have approximately flat circular velocity curves, spheroids tend to be denser and to have declining Vc curves. 



We show the ej distribution in the right-hand panels of Fig.[T] 
for three simulated galaxies, chosen to illustrate three representa- 
tive cases. The top panel corresponds to a galaxy where most stars 
are in coplanar, nearly circular orbits, hence the sharply-peaked dis- 
tribution near = 1. The bottom panel corresponds to a spheroidal 
galaxy where ordered rotation plays little role; the ej distribution is 
broad and centered around zero. The middle panel corresponds to 
an intermediate case, where a non-rotating bulge of stars is sur- 
rounded by a well-defined thin disk. A simple quantitative mea- 
sure of morphology can therefore be constructed by the fraction 



of stars with circularities exceeding a fixed fiducial value, such as 
/fe > 0.5). 

Although conceptually simple, ej distributions are not easy to 
relate to observation, so a simpler quantitative measure of morphol- 
ogy is desirable. One alternative is the fraction of kinetic energy 
invested in ordered rotation, 

trot ~ 1 for disks with perfect circular motions, and is 1 for 
non-rotating systems. As Fig. [T] makes clear, Krot correlates ex- 
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Figure 4. The kinematic moiphology parameter, Krot , versus a number of parameters characterizing the properties and assembly history of each galaxy and 
its surrounding halo. On the left, from top to bottom, t5o% is the half-mass halo formation time, in Gyrs; AA/i„j„j is the maximum fraction of the final halo 
mass assembled in the single largest merger event after 2 = 3; and A' is the dimensionless rotation parameter (eq. |2)- On the right, the galaxy formation 
"efficiency", '(ygai,* = A/gai/(/bar ^^^^200); /acc is the fraction of accreted stars (i.e., stars formed in galaxy progenitors other than the main one) and /i^ot is 
the fraction of stars bom out of gas that went through the "hot phase" (i.e., Tmax > 10^'^ K). Correlation coefficients for each panel are given in Table|2] See 
text for further details. 



tremely well with the fraction of stars with e > 0.5. In what fol- 
lows, we shall use Krot to rank galaxies according to the importance 
of their rotationally-supported components. 

For convenience, we shall hereafter refer to galaxies with 
Krot < 0.5 and n^ot > 0.7 as spheroid- or disk-dominated, re- 
spectively. The first group makes up ~ 25% of the sample; the 
second group comprises another ~ 30%. The remainder consist of 



intermediate types where both rotation and velocity dispersion play 
a comparable structural role. It is important to note that our sample 
contains galaxies spanning a wide range in Krot, from "pure disk" 
systems with a negligible fraction of stars in counter-rotating orbits 
(i.e., < 0) to spheroids with little trace of rotational support. 
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3.2 Examples of galaxy morphologies 

Fig. |2] shows four examples chosen to illustrate the structure of 
galaxies with various values of Krot • The panels show edge-on and 
face-on projections of each galaxy, colored by stellar surface mass 
density on a logarithmic scale. Fig. [3] shows circular velocity pro- 
files and (face-on) stellar surface density profiles. The degree of ro- 
tational support increases from left to right: the leftmost and right- 
most are spheroid- and disk-dominated systems, respectively, while 
the two middle ones are intermediate-type objects. Labels in each 
panel indicate, for each galaxy, the stellar mass within the radius, 
''gal = 0.15 r2oo, used to define the central galaxy. Table [T] lists 
some physical parameters of galaxies A-D. 

Figs. |2] and [3] together with Table [T] show that simulated 
galaxies have several properties in common with nearby ellipti- 
cals and spirals. Spheroid-dominated galaxies are gas-poor, dense 
stellar systems with declining circular velocity curves, whereas 
disk-dominated galaxies are richer in gas, more spatially extended, 
and have nearly-flat circular velocity curves. Interestingly, from 
the point of view of surface density profiles, spheroids are single- 
component systems well approximated by de Vaucouleurs' B}^'^ 
law (dashed magenta lines in the bottom panels of Fig.O. Disk- 
dominated systems, on the other hand, have more complex profiles, 
with a central R^^'^ spheroid surrounded by an exponential compo- 
nent that increases in importance in step with Krot ■ Like most spi- 
rals, they are well approximated by the sum of a de Vaucouleurs' 
spheroid (dashed magenta lines) and an exponential law (dashed 
cyan lines). 

As discussed in detail by McCarthy et al (in preparation), these 
similarities with observation actually extend to quantitative com- 
parisons with observed scaling laws, such as the Tully-Fisher rela- 
tion or the Fundamental Plane. The agreement between simulated 
galaxies and observation is encouraging, and suggests that the ori- 
gin of the morphological diversity of simulated galaxies can pro- 
vide insight into what determines the relative importance of disks 
and spheroids in real galaxies. 

A few caveats should also be mentioned. A wide variety of 
morphological types does not automatically guarantee reproduc- 
ing the right morphological mix of galaxies in this mass range. 
Indeed, a casual inspection suggests that spheroids might be over- 
represented in our sample. However, because most morphological 
classifications are based on photometric data, exploring this issue 
in detail would require synthesizing "observations" of the simu- 
lated galaxies in various bands and analyzing them in the same 
way as observed samples, which is bey o nd the scope of this work. 
Moreover, as discussed in ICrain et all ( |2009|) . the GIMIC galaxy 
stellar mass function differs substantially from the observed one, 
implying that such exercise would be inconclusive, regardless of 
its results. Our main goal is thus not to test the viability of the par- 
ticular star formation/feedback implementation adopted in GIMIC, 
but rather to learn about the various mechanisms responsible for 
the origin and relative importance of rotationally-supported versus 
dispersion-supported stellar components in a galaxy. 
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Figure 5. The fraction of stars bom from gas that went through the "hot 
phase" (i.e., Tmax > lO^-S K) versus the median formation time of stars in 
the galaxy, tgpo^ ■ Symbols encircled in grey correspond to systems where 
more than 10% of stars were accreted in merger events. Each galaxy is color 
coded according to its morphology: red dots, green asterisks and blue trian- 
gles correspond to spheroid-dominated, intermediate, and disk-dominated 
systems, respectively, as labeled by the Krot parameter Histograms show 
the distribution of median formation time for each of these three groups. 
The tight correlation between /jiot and *5o% suggest that gas able to reach 
the "hot phase" take longer to accrete and to be transformed into stars that 
gas accreted "cold". Late gas accretion clearly favors the assembly of stellar 
disks. 



4 THE ORIGIN OF SIMULATED GALAXY 
MORPHOLOGIES 

4.1 Dependence on dark halo properties 

As discussed in Sec.[T] stellar disks are expected to form at the cen- 
ters of halos with quiet recent accretion histories and high angular 
momentum. Halos that have been relatively undisturbed by recent 
major mergers tend to form earlier, so we may also expect stellar 
disks to inhabit halos with early formation times. 

We explore this in Fig. |4] where the left panels show the de- 
pendence of Krot on (i) the halo formation time, t5o% (when the 
most massive halo progenitor reaches half the final halo mass); on 
(ii) the fraction of halo mass accreted in the single largest merger 
after z = 3, AAfimm, and on (iii) the dimensionless rotation pa- 
rameter. 



^2 iW'2ooV2oor2oo ' 

where J is the total angular momentum of the halo jBullock et al.l 
I2OOII) . 

None of these parameters correlates strongly with galaxy mor- 
phology (see Table [2J. Disks form in halos with low and high spin 
parameter; in halos that collapse early and late, and even in halos 
that have accreted a substantial amount of mass in merger events. 
The same applies to spheroids, except perhaps for a weak tendency 
to prefer halos with slightly lower- than-average A'. 

Fig. |4] also shows that major mergers are uncommon during 
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black; in-aitu stars red: accreted stars 




Figure 6. Fomiation time of stars (expressed in terms of redshift) versus 
the circularity parameter, ej , measured at z = 0, for the four galaxies illus- 
trated in Fig. [2] Stars formed in the main progenitor (in situ) are shown in 
black, accreted stars in red. The magenta curve tracks the median circularity 
as a function of formation time. The two spheroid-dominated galaxies show 
signs of episodic star formation events that lead to the presence, at 2 = 0, 
of stellar populations with distinct angular momentum properties. Disks, on 
the other hand, tend to form more gradually over time and to be dominated 
by a single population with coherently-aligned angular momentum. 



the formation of halos in the narrow mass range considered here; 
0.5 < M2oo/10^^ Mq < 1.5. Most systems have accreted 
less than 20% of their final mass in a single merger since z = 3, 
and these events seem unrelated to the morphology of the central 
galaxy at 2 = 0. 

Finally, morphology also seems unrelated to the fraction of 
baryons within the virial radius that collects to form the galaxy. 
This is illustrated in the top-right panel of Fig. [4] where we plot 
/trot vs the galaxy formation "efficiency" parameter, r?gai,* = 
Mgai,,/{fh3.i M200), where /bar = ^b/^m = 0.175 is the uni- 
versal baryon fraction. Although we consider halos in a narrow 
mass range, the efficiency of galaxy formation varies from 20% 
to 70% (with an average of (??gai,*) = 40%) and appears to have 
little influence on the morphology of the central galaxy although 
there is a weak tendency for disky objects to prefer lower values of 

Thus, contrary to simple expectations, spheroids can form in 
quiescent halos, and disks can form in halos of scant angular mo- 
mentum content. Simple predictions of the morphology of a galaxy 
based on the properties and assembly history of its surrounding 
dark halo will thus often be wrong. 



^ We note that abundance matching models (see, e.g.. lMoster et alJboiOl ; 
iGuoetalJlIoToh demand an even lower galaxy formation efficiency to 
match the galaxy stellar mass function, which would weaken even further 
the link between the properties of a central galaxy and that of its surround- 
ing halo. 



4.2 Dependence on galaxy history 

Galaxy mergers can still in principle play a role in determining mor- 
phology, if their importance is underestimated by the halo merger 
parameter AMimm. Indeed, galaxies take longer to merge than ha- 
los do, and, due to the large scatter in galaxy formation efficiency, 
the mass ratio of galaxy mergers may differ substantially from that 
of their surrounding halos. 

We examine the importance of accretion on morphology more 
explicitly in the middle-right panel of Fig.|4l where we plot Krot vs 
/acc, the fraction of stars accreted by the galaxy; i.e., those formed 
in systems other than the main progenitor of the galaxy. This is a 
direct measure of the importance of accretion events in the build- 
up of the galaxy. Two points are worth noting here: most galaxies 
form the majority (> 90%) of their stars in-situ, and there is no 
correlation between Krot and /acc- The accreted fraction exceeds 
25% in only 5 of our 100 simulated galaxies; overall, accretion 
events just seem to bring in too few stars to play a significant role 
in the morphology of our simulated galaxies. 

An interesting clue is provided by the thermodynamic history 
of the gas before it is transformed into stars. This may be estimated 
simply by tracking every stellar particle back in time and by record- 
ing the maximum temperature, Tmax, reached before the particle 
accretes into the galaxy and becomes eligible for star formation. 
If Tmax exceeds 10^'^ K, then in all likelihood it was accreted by 
gradual c ooling f rom a shock-heated, nearly hydrostatic gas corona 
(see, e.g.. lCrairi e t al. 2010; van de Voort et al. 201 1). 

The fraction of stars, /hot, whose parent gas particles went 
through this phase correlates well with /trot, indicating that accre- 
tion of gas from the "hot phase" favours the formation of disks 
(see bottom-right panel of Fig. |4}. No disk-dominated galaxy (i.e. 
/trot > 0.7) forms unless /hot exceeds 40%. This is intriguing, 
since it runs counter recent proposals that "cold flows", i.e., gas 
that gets accreted directly into the galaxy without going through 
the hot phase, might promote the for mation of exte nded disks (e.g., 
iKeres et al J 120051 ; iDekel et al.l 120091 ; ISrooks et al J 12009). If any- 
thing, our simulations suggest the opposite; the majority of stars 
in spheroid-dominated galaxies originate in gas that accretes cold. 

One reason for this result is illustrated in Fig. [5] which shows 
that there is a tight correlation between the fraction of stars bom 
from gas that cooled from a hot corona and the median formation 
time of stars in a galaxy, t5o%- Heating gas to a hot corona be- 
fore cooling delays its accretion and favors the late assembly of a 
galaxy: the larger /hot the later stars form. Enhancing recent star 
formation promotes the formation of disks and facilitates their sur- 
vival until the present. This scenario, although appealing, seems 
incomplete, given the presence of several spheroid-dominated sys- 
tems that form despite high /hot, late t5o%, and in the absence of 
significant merger activity. 

4.3 Dependence on spin alignment 

Our results so far suggest that the morphology of GIMIC galax- 
ies is linked largely to internal mechanisms operating in individ- 
ual galaxies rather than to accretion-driven transformations. This 
has been anticipated by semi-analytic models of galaxy formation, 
where secular evolution driven by "disk instabilities", is thought 
to be an important formation path for spheroids. These instabil- 
ities are assumed to be triggered when the self-gravity of a disk 
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Figure 7. Projected particle distribution neai' turnaround time, z = 3.5, of baryons that collapse to form, at 2 = 0, galaxies B and D shown in Fig.|2] Stars 
already formed are shown in red, particles still in gaseous form in blue. Box sizes are in physical units. Concentric circles enclose 20%, 50%, and 95% of the 
mass, and arrows indicate the angular momentum of all material enclosed within each radius. Arrow lengths are normaUzed to the total value, which defines 
the z axis of the projection. Each panel is normalized separately, so that jo. 95 has equal length in both. Note the misalignment of the angular momentum of 
various parts of the system for the spheroid-dominated galaxy B. Angular' momentum is more coherently acquired in the case of the disk-dominated galaxy D. 




Figure 8. Angle between the angular momentum enclosed within a given mass fraction, m / rrttot , and the total spin of the system measured at the time of 
maximum expansion (turnaround) for twenty galaxies in our sample: the ten systems with highest and lowest values of Krot at 2 = and /acc < 0.1, 
respectively. The restriction in /acc is included in order to focus on systems unaffected by merger events. By definition, all curves approach unity as m — > 
mtot. Panels on the left correspond to all baryons within the galaxy radius, Tgai, at 2 = 0; those on the right to dark matter halo particles that are within r200 
at 2 = 0. Note the strong misalignments between different parts of the system for galaxies that are spheroid-dominated at 2 = 0, and the smooth alignment 
characteristic of the turnaround configuration of systems destined to become disks. 
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Figure 9. Same as left-hand panel of Fig. [8] but averaged over all galaxies 
grouped in bins of Krot (at 2 = 0), as labeled. Note that the more disk- 
dominated a galaxy is at present, the more coherently aligned the spin is 
at the time of turnaround. The inset panel shows (cos(9)} for individual 
systems, averaged over all mass shells at turnaround, as a function of Krot ■ 
Galaxies where accretion events have played a minor role, i.e., /acc < 0.1, 
are shown with solid circles, the rest with crosses. 



exceeds a particular threshold. For fixed halo mass, as in our sam- 
ple, this should lead to noticeable correlations between the mass 
of the galaxy and the importance of the spheroid. However, as we 
discussed above, no such correlation is apparent. 

Recalling the relation between /hot, ^50%' and ''^rot, we look 
at the star formation history for further clues. Fig. |6] plots the for- 
mation redshift, Zform, of stars in the four galaxies shown in Fig.|2] 
versus the "circularity" parameter, ej, measured at z = 0. Points 
in black correspond to stars formed in-situ (i.e., within the main 
progenitor) while those in red are accreted stars. 

This figure shows that, as expected, star formation proceeds 
gradually in disk-dominated systems. In spheroids, however, stars 
form in separate episodes that leave behind stellar "populations" 
of different ages and distinct angular momenta. Since most stars 
form in-situ, these populations in galaxies A and B are likely caused 
by gas accretion events where the net angular momentum of one 
event is misaligned with the others. On the other hand, disks tend to 
form out of accreted gas that shares a common angular momentum 
direction. 



These considerations suggest that the final morphology of a 
galaxy is imprinted early on, since the spin of the material destined 
to form a galaxy is acquired at the time of maximu m expansion 
and changes little in the absence of merging (see, e.g.. lWhitel 19841 : 
iNavarro et alj2004l) . We illustrate this by studying the angular mo- 
mentum of galaxies B and D at z = 3.5, which roughly corre- 
sponds to the time of turnaround of both systems. Fig. [7] shows 
the spatial distribution of all baryons that will end up within rgai 
at z = 0. Net angular momentum is acquired through the inter- 
play between the inertia tensor of the mass distribution and the 



ll996VPorc iani et al.l20023lbl:lNavarro et alj2004) . For highly non- 
uniform spatial distributions, where the principal axes of the inertia 
tensor can change direction abruptly, this effect may cause the net 
angular momentum of different parts of the system to flip and mis- 
align. 

We see from Fig. [7] that this is indeed the case for galaxy B. 
Here the arrows indicate the direction and magnitude, at 2 = 3.5, 
of the specific angular momentum of the inner 20%, 50% and 95% 
of the baryons that end up in the galaxy al z — 0. The length of the 
arrows is normalized to the total, which is chosen by construction 
to coincide with the z axis of the projection. The angular momenta 
of different parts of the system are clearly misaligned: the angle 
between jo. 2 and jo. 5 is 85 degrees, and that between jo. 5 and jo. 95 
is ~ 100 degrees. Since gas further out in Fig. [7] takes longer to 
accrete, newly assembled material will be misaligned with the rest, 
leading to the formation of distinct populations of stars (shown in 
Fig.|6ll that will tend to destabilize any existing disk and to cancel 
out the net angular momentum of the system; leaving in place a 
slowly-rotating stellar spheroid. 

On the other hand, the spins of different parts of the system 
are very well aligned in the case of the disk-dominated galaxy D, 
as shown in Fig. |7] This coherence in the angular momentum at 
turnaround allows newly accreted material to settle into a stable 
disk where star formation can proceed gradually and smoothly. 

We show in Fig. [8] that this result applies to the majority of 
spheroid- and disk-dominated galaxies in our sample. Here we plot, 
for the ten galaxies with highest and lowest values of Krot where 
accretion has played a minor role (/acc < 0.1), the cosine of the 
angle 9 between the angular momentum of a given enclosed mass 
fraction m/mtot and the direction of the total spin of the system. 
By construction, each curve in Fig. [8] is constrained to approach 
unity as the enclosed mass approaches mtot. 

It is clear from this figure that different regions of systems 
destined to form spheroid-dominated galaxies have, at turnaround, 
large misalignments in their acquired spin. Indeed, in many cases 
the inner regions counterrotate (i.e., cos(6) < 0) relative to the 
outer regions of the system. This is not the case for systems that 
become disk-dominated which, in general, show coherence in the 
alignment of the spin axis. In these cases, the enclosed specific 
angular momentum increases roughly linearly with enclosed mass 
fraction. As Fig. [9] shows, the same result applies to all galaxies 
in our sample: despite the large scatter, on average, the degree of 
alignment at turnaround increases gradually with the importance of 
the disk component in the morphology of a galaxy at z — 0. 

The right panels in Fig. [8] show that a similar assessment ap- 
plies to the dark matter halo surrounding these galaxies. The halos 
of spheroids also show, at turnaround, stronger misalignments than 
the halos that host disk galaxies at z = 0. This is encouraging, 
since it implies that it might be possible to use the angular momen- 
tum properties of a dark matter halo at turnaround to "predict" the 
morphology of its central galaxy at z = 0. We emphasize, how- 
ever, that the trends we highlight here, although well defined, are 
relatively weak, so the correspondence between early halo prop- 
erties and final galaxy morphology is likely to apply statistically 
rather than to individual systems. 

Quantitatively, the dependence of present day morphology on 
spin alignment at turnaround is shown in the inset panel of Fig. [9] 
Here we show, for each individual system, Krot versus the average 
cosine of the angle, at turnaround, between the angular momen- 
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Figure 10. Fraction of stars born from gas accreted from the "hot phase" as 
a function of the ahgnment parameter (cos(9)) at the time of turnaround. 
Color-coding is the same as in Fig. [5] Note that, in general, disk-dominated, 
spheroid-dominated, and intermediate systems occupy different regions in 
this diagram. See text for further discussion. 



turn enclosed by different mass shells and the total. The trend is 
clear: systems with better-aligned spins at turnaround tend to be 
more disk-dominated at present. The trend is even stronger when 
considering only systems where mergers have played a minor role 
(.face < 0.1, solid points). The correlation coefficient is = 0.49 
(with significance A --^ 2 x 10~®) and confirms our earlier con- 
clusion that significant misalignment in the distribution of angular 
momentum can lead to the formation of spheroid-dominated sys- 
tems in the absence of significant merger events. 



4.4 Spin alignment vs mode of accretion 

The results of the above subsections suggest that there are (at least) 
two mechanisms responsible for the morphology of GIMIC galax- 
ies: the alignment of the spin acquired by various parts of a galaxy 
and the mode of gas accretion. Fig. [TO] shows that these two effects 
are approximately independent of each other: there is no obvious 
correlation between the fraction of gas accreted from the hot phase 
and the alignment parameter (cos(^)). However, disk-dominated 
systems (blue triangles) separate from spheroid-dominated ones 
(red circles) fairly neatly in this plane, suggesting that Krot is de- 
termined by a combination of the two mechanisms. 

Large /hot favors disk formation: gas shock-heated into a 
nearly-hydrostatic corona of hot gas is forced to homogenize its ro- 
tational properties before accretion, providing the forming galaxy 
with a gradu al supply of gas that shares the same spin axis (see, e.g., 
iBrook et a l. 2012). Disk-dominated systems then form when /iiot 
is large, unless the spin misalignments are large enough to disrupt 
them: no disk-dominated system forms when (cos(S)) < 0, even 
for relatively large values of /iiot • 

Spheroids, on the other hand, form primarily when cold gas 



accretion prevails (/iiot < 0.4): gas that flows along distinct fila- 
ments cannot interact hydrodynamically before accretion and will 
often have misaligned net spins. Each accretion event then results 
in the formation of a "population" of misaligned stars that will tend 
to destabilize any existing disk and to cancel out the net angular 
momentum of the system, leaving in place a slowly-rotating stellar 
spheroid. Intermediate systems, in general, result when the delete- 
rious effects on stellar disks of cold accretion events are mitigated 
by well-aligned spins. 

In other words, spheroid-dominated galaxies in our sample 
do not originate from disk instabilities triggered by self-gravity, 
as envisioned by semi-analytic models, but rather by the accre- 
tion of gas that settles on off-axis orbits relative to ear lier accreted 
mater i al. This has been seen in earlier work (e.g.. iBrook et al.l 
l2008l : IScannapieco et al.ll2009l) . and might be related to sudden 
c hanges in the orienta tion of the dark matter halos as discussed 
m lBett&FrenkH201lh . Its relevance to the formation of the whole 
class of spheroidal galaxies, in the absence of merging, has not yet 
been recognized and emphasized. 

We hasten to add that the importance of this mechanism for the 
formation of spheroids might depend on halo mass, and that we ex- 
plore only a narrow range here: 0.5 < Af2oo/10^^ Mq < 1.5. 
This caveat might be particularly relevant in the case of the most 
massive s pheroids, where merging likely plays a more important 
role (e.g. IParrv et al.ll2009l : lOe Lucia et aUl201ll : iFeldmann etaTI 
I2OIOI) . Caution must also be exercised when extrapolating the link 
between morphology and the fraction of stars born from hot accre- 
tion. Disk-dominated galaxies might still form out of cold accretion 
if, for example, the most recent episode of accretion supplies most 
of the mass of the system. Repeated cold accretion events may hin- 
der disk formation, but a single major event may very well facilitate 
it. 



5 SUMMARY 

We use gasdynamical cosmological simulations of galaxy forma- 
tion to study the origin of different galaxy morphologies in the 
ACDM cosmogony at redshift z = 0. The GIMIC simulation series 
covers a large volume and has a resolution high enough to study 
the structure and kinematics of the stellar components of 100 cen- 
tral galaxies in Milky Way-sized halos. We focus our analysis on 
the origin of galaxy morphology, somewhat narrowly defined as 
the relative importance of rotational support vs velocity-dispersion 
support (the disk-to-spheroid ratio) in the structure of the galaxy. 
Our main results may be summarized as follows: 

• The simulated galaxies span a wide range of morphological 
types, from rotation-free spheroids to almost pure disk galaxies 
where fewer than 5% of all stars are in counterrotating orbits. Disks 
have roughly exponential stellar surface density profiles and flat 
rotation curves, whereas spheroids are dense stellar systems that 
can be approximated by de Vaucouleurs' R^^'^ profiles. The re- 
semblance with real galaxies suggests that it should be possible to 
gain insight into the origin of galaxy morphology by studying the 
mechanisms responsible for the relative importance of disks and 
spheroids in GIMIC galaxies. 

• The morphology of simulated galaxies seems mostly unrelated 
to the spin or assembly history of their surrounding dark matter 
halos. Most stars form in- situ and comprise on average about ~ 
40% of all available baryons in the halo. Most baryons in a halo 
therefore end up not making part of the central galaxy, which helps 
to explain the weak correlation between the properties of halos and 
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those of central galaxies. Contrary to simple expectations, disks 
form in halos with low and high angular momenta, and spheroids 
form even in galaxies where most stars form in-situ, suggesting a 
formation path for spheroids that does not rely on merging. 

• The star formation history provides an interesting clue to the 
origin of morphology. Disks tend to have young stars, and to form 
gradually over long periods of time. This is because gradual cooling 
from a hot corona delays the accretion of gas and promotes late star 
formation. 

• Star formation in spheroids proceeds episodically, leaving be- 
hind populations of stars of similar age but distinct kinematics. 
These populations originate from the accretion of gas whose an- 
gular momentum is misaligned relative to that of earlier-accreted 
material. The misalignment destabilizes any pre-existing disk, 
prompts the rapid transformation of gas into stars, and reduces the 
net rotational support of the system. 

• Since angular momentum is largely acquired at the time of 
maximum expansion of the material destined to form a galaxy, a 
good indicator of morphology at ^ = is the coherence in the 
alignment of the net spin of various parts of the system at the time 
of turnaround. Spheroid-dominated galaxies form in systems where 
misalignments are substantial whereas disks form in systems where 
the angular momentum of all mass shells is roughly aligned. 

• The final morphology of a galaxy results from the combined 
effects of spin alignment and of hot/cold gas accretion. Disk- 
dominated objects are made of stars formed predominantly in situ, 
and avoid systems where most baryons were accreted cold, or those 
where spin misalignments are extreme. On the other hand, direct 
filamentary accretion of cold gas, especially when accompanied by 
substantial spin misalignments, favours the formation of slowly- 
rotating spheroids, which may thus form even in the absence of 
mergers. 

Our results suggest a new scenario for the origin of ~ i* stel- 
lar spheroids that does not rely on merging. This scenario, once 
developed more thoroughly, should offer a number of predictions 
falsifiable by observation. For example, the episodic nature of star 
formation in spheroids envisioned here is expected to leave be- 
hind overlapping populations of stars of distinct age, kinematics, 
and, possibly, metallicity) that survive to the present because of the 
paucity of mergers. We plan to explore the observational signatures 
of these populations in future work. 

The scenario we propose here also offers clues to the origin of 
pure disk galaxies. A number of our simulated galaxies have vir- 
tually no "classical" spheroid, with fewer than 5% of their stars in 
counterrotating orbits. These galaxies form either in systems where 
spin alignment is extraordinarily coherent or where most of the 
baryons in the galaxy get accreted late from a hot corona. Although 
at this time limited numerical resolution precludes a more detailed 
study, we plan to use these clues to resimulate some of these sys- 
tems at higher resolution with the goal of shedding light into the 
origin of bulgeless galaxies. 

We emphasize that, although coherent spin alignment at early 
times is clearly an important clue, it should be considered as one 
ingredient of the complex process that determines the morphol- 
ogy of a galaxy. Strong feedback, for example, may expel baryons 
from galaxies and cycle them through a hot corona before they 
get re-accreted and turned into stars, potentially erasing the spin 
alignment dependence at turnaround we report here. Furthermore, 
aligned spins in the accreting gas might not be enough to ensure 
the survival of a stellar disk, especially if the dark matter halo is 
strongly triaxial and its principal axes are not coincident with the 



Table 1. Summary of main properties for Gal A-D in Fig. |2] and [3] Rows 
correspond to the virial mass Af20o; galactic mass in stars A/str; gas Afgas 
and gas fraction /gas; peak circular velocity Vlnax; the circular velocity 
measured at the galactic radius Vc{r = fgai); the degree of rotational sup- 
port Ki ot ; and the fraction of the stars in counter-rotating orbits fc . 



Property 


Gal A 


GalB 


GalC 


GalD 


Af20o[10i2h-iAfQ] 


1.18 


0.98 


0.77 


1.04 


Mstr(r < rgai)[10iO/i-iAfQ] 


8.98 


7.25 


5.42 


6.31 


Mgae(r < rgai)[lOiO/i-lAf0] 


0.51 


0.72 


1.00 


1.86 


/gasCAfgas/CAfgas + Afstr)) 


0.05 


0.09 


0.15 


0.23 


Kiax [km/s] 


471 


383 


300 


280 


Vc(r = rgai) [km/s] 


187 


172 


160 


178 


^rot 


0.31 


0.47 


0.61 


0.76 


/c 


0.48 


0.33 


0.12 


0.06 



Table 2. Spearman rank correlation coefficients, r^, between ftrot and the 
halo/galaxy properties shown in Figs.|4]and|9] The second column shows 
the two-sided significance of its deviation from zero, as computed by the 
IDL subroutine r-correlate: smaller values indicate more significant 
cortelations. 



Property 










A 


*50% 


0.05 






0. 


60 




-0.06 






0. 


57 


A' 


0.29 


2.S 


! X 


10- 


-3 


Vgal,* 


-0.30 


1.9 X 


10^ 


-3 


face 


0.06 






0. 


55 


/hot 


0.55 


i.e 


; X 


10- 


-9 


< cos (9) > 


0.41 


2.0 X 


10 


-5 



disk. Finally, although mergers are rare in the mass range we ex- 
plore here, they likely play a more important role in the forma- 
tion of more massive spheroids. Until simulations can reproduce 
not only the properties of individual systems, but the full statistical 
distribution of galaxy morphologies and their dependence on mass 
and environment, it is likely that a full understanding of the origin 
of galaxy morphology will remain beyond reach. 
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